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— 1  This  paper  discusses  autoregressive  random  field  (ARF) 
models  and  derives  a  unilateral  representation  driven  by 
uncorrelated  noise. 
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1.  Introduction 


In  this  paper  we  shall  deal  with  a  subset  of  stationary 
(wide  sense)  processes  with  absolutely  continuous  spectral 
distributions  which  are  rational  functions  of  the  two  quan- 
tities  e  ,e  .  More  precisely  we  shall  study  the  process 


X[m  njfeR  ,  [m(n]6ExB  on  an  infinite  lattice,  with  covariance 


structure 

E  (X 


[m+s,n+t]X  [m,n]^ 


i  „  -ise.-ite 

-±2 f  f  e  f(8 ,,0 2)dQ1dQ2, 

4tt  —it— it 


(1.1) 


and  zero  mean. 


-1 


We  assume  exists  and  is  finite  at  every  (0^,©2), 


and 


*  f(8  ,©  )  1  =  (a ,Q  I  a ,  .cos (m, 0 , +m_6_) ) 

1  2  10,01  ^l'^1  2  2 


(1.2) 


a[0,0],a[mi,m2]  are  p  x  p  matrices  satisfying  • 

V*  is  the  complex  conjugate  transpose  of  the  vector  V.  np 
denotes  the  deleted  p  x  p  neighborhood  of  10,0],  that  is, 

{  [m1,m2]  :  (mj^l  sp,  |m2|  s.p,  (m1,m2]^[0,0]  }. 

Models  of  this  type  have  been  used  as  models  of  texture 
images  [1,2].  in  the  case  where  Xr  ,  is  a  Gaussian  process, 
it  can  be  shown  [3]  that  Xj  j  is  a  Gauss- Markov  process  with 
respect  to  NP;  that  is, 


(1.3) 


P‘ob(X[m,n]  'X[s,t]'  = 

Prob(X[m,n]  |x[s,t]  '  U,t]€[m,n]+Np) 

In  fact,  the  process  with  spectral  density  f(0]/e2)  satisfies 
the  conditional  model 

aIn,1.m2]x[m-ml,n-m2]»  <1-4' 


Conditional  models  of  this  type  have  been  found  useful  in 
the  modeling  of  spatial  patterns  [7].  It  is  also  known  (see, 
for  example,  page  26  of  [7])  that  no  finite  one-sided  repre¬ 
sentation  for  this  model  exists  of  the  type  (with  S  finite 
subset  of  2Z  x  22) 


blO,0]X(m,n]  +  [ni^ro2]€Sb[m1,m2)X[m-m1,n-m2]=Z[m,n] 

where  Z nj  is  a  process  of  uncorrelated  noise. 

The  purpose  of  this  paper  is  to  show  that  the  collection 
of  spectral  representations  of  the  process  Xr  ,  along  one 

l  •  9  •  i 

of  the  coordinates  is  representable  as  a  one-sided  finite 


order  "time  series"  model  along  the  other  coordinate.  Thus, 
in  this  sense  it  is  seen  that  all  ARF's  have  a  "causal" 
representation . 

This  method  of  producing  a  one-sided  representation  can  be 
contrasted  with  the  so-called  NSHP  (non- symmetric  half  plane) 
representation  of  [3]  and  [6]. 


2 .  A  Unilateral  Representation 


We  consider  the  process  Xr  .  with  spectral  density 

L  •  r  •  J 

f  (9,  ,9,)  =  [afn  n1  +  E  a,  ,008(111,0+111-0-)] 

1  2  10 »0]  r_  _  -.,„P  [m,  ,nul  11  2  2 


-1 


10'01  'V"2I 


(2.1) 


which  is  a  p-th  order  autoregressive  process. 

-i©! 

Let  z=e  ,w=e  ,  and  rewrite  the  above  equality  as 
f (01 , 02) ~^“ao +ai M  z+* • • +*p (w) zp 

+a* (w) z~1+. . ,+a*  (w)z-p. 

For  each  fixed  w  we  can  consider  f(0^,02)  as  a  spectral 
density  in  0^.  We  next  produce  a  causal  factorization  of 

^°1'°2^  ;*’n  t*ie  ^orm 

f(01,02)’1=  (2.2) 

^  » 1  It  **13  n 

[Cq  (w)  +C^  (w)  Z  +...+Cp(w)Z  P]  [Cq  (wJ+Cj^  (w)  Z+.  .  ,+Cp(w)  ZF]  , 


”i02 

where,  for  each  w=e  ,  c0(w)+c.(w)£+...+c  (w)£p  has  no 

roots  inside  and  on  the  complex  unit  circle  |£|=1.  ([4], 

page  65) . 

We  next  consider  the  spectral  representation  of  the 
process  Xf  .  along  the  second  coordinate: 

l  •  t  •  i 


Xln.ml*.^i,n6dyn(e>' 


(2.3) 


where  Yn ( 0 )  is  the  spectral  representation  of  the  process 

X[n,.].  (t5J'  Pa^e  481). 

Next  expand  each  of  cQ (w) , . . . ,cp (w)  in  a  Fourier  expansion 


I 

l 


c  .  (w)  =?  e 
3  k=-o° 


ik  6 


2 


Then  (14],  page  61)  the  process 


satisfies  the  autoregression 


p  00  A 

j=o  k=-»Ij/klXln"j#m+lcl=Z[n,ml 


where  Z 


(2.4) 

is  an  uncorrelated  white  noise  process.  Let 


Wn(0)  be  the  spectral  representation  of  the  process  z[n#<j: 

7  f11  imO,t.  (2.5) 

[ n,m]  -ir  n' 

We  conclude  with  the  following 

Theorem:  Let  <Yn(0)},  (Wn(0)}  be  the  spectral  representations 
of  the  processes  defined  above.  They  satisfy  the  stochastic 


differential  equation 


P  — i0 

l  ck(e  1O)dYn_k(0)=dWn(0)  (2.6) 

k-0 


Proof:  In  the  above  autoregressive  representation  we  sub¬ 
stitute  the  spectral  integrals  and  get  (after  combining  terms) 


Vm:  /"{I  l  cr.  el(m+k)0dY  . (6) -eim0dw  (0) ]  *0 

-ir  j=0  k=-~  13  ,KJ  n"3  n 


Factoring  out  eim0  we  have 

7T  p 

Vm:  /  eim0U  c . (e~10 ) dY„  . (0)-dW  (0)}  =0 
-ir  j=0  3  n-3  n 

As  any  continuous  function  f  (0)  ,0fc[-Ti,Tr)can  be  approximated  in 

mean  square  by  linear  combinations  of  eim0 


,  the  result  follows 


3.  The  finite  version. 


The  above  calculation  can  be  carried  out  in  the  case 
where  we  have  a  finite  number  of  values 

X[n,0] - - X[n,M-l] 

in  the  vertical  direction. 

Let  iji  =e^Tt3m^.  The  finite  versions  of  the  above  spectral 
representations  are  as  follows: 


M-l 

X  r  ,  =  E  i b  km„ ,  ,  > 

[i,m]  k=0  yM  V(n,k) 


or 


1  M-l 


-km 


4Y(n,k)— 


[n,  j] 


(3.1) 

(3.2) 


Tnat  is,  AY (n, • )  is  the  finite  Fourier  transform  of  the  data 

Xr  ,.  Similarly 
in,  •  j 


1  M-l 

Aw (n,k)  = -  E  w.-km  _ 

M  j=0  VM  Z  [n, j ] ’ 


(3.3) 


The  finite  analogue  of  the  above  Theorem  is 
M-l  ink  p  -k 

1  {  L  b.  (ipM  )  AY  (n- j  ,k)  -  AW(n,k)}  =0 

k=0  M  j=0  J 

concluding  with 

p  -k 

E  b.  (ij<MK)AY(n-j,k)=AW(n,k) . 
j=0  3  M 


(3.4) 
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